More on Stochastic and Variational 
Approach to the Lax-Priedrichs Scheme 
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The author presented a stochastic and variational aspect of the Lax-Friedrichs 
scheme applied to hyperbolic scalar conservation laws [18| . We will extend the 
results of [18] on the Lax-Friedrichs scheme, showing its time global stability; its 
long time behaviors; its error estimates. The proofs essentially rely on calculus 
of variations in the Lax-Friedrichs scheme and viscosity solutions of Hamilton- 
C$ . Jacobi equations corresponding to the hyperbolic scalar conservation laws. We 

also provide basic facts that are useful in numerical analysis and simulation of 
the weak KAM theory. As an application, we show rigorous treatment of finite 
difference approximation to KAM tori. 
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S-J '■ 1 Introduction 

> : 

We investigate the Lax-Friedrichs scheme applied to initial value problems of hyperbolic 
H ! scalar conservation laws with a constant c 

(1.1) u t + H{x,t,c + u) x = 0. 

There is the huge literature on stability and convergence of the scheme. The standard 
technique is based on the L 1 -framework with a priori estimates and the compactness 
of functions of bounded variation, where mesh size independent boundedness of both 
difference solutions and their total variation must be verified. Since the Lax-Friedrichs 
scheme is very simple, one can successfully analyze details of approximation, particu- 
larly in the case where a flux function is of the simple form H(x,t,p) = H(p). We 
refer to [5], [19], [15] and the works cited there. However, in the case of a general 
flux function depending on both x and t, the problem becomes much harder and often 
requires unpleasant assumptions. The results of the general case first appear in [T4] . 
where stability and /^-convergence are proved with a restricted time interval that is 
determined by the growth of H(x,t,p) with respect to p. In [13], time-global stability 
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and L 1 -convergence within arbitrary time intervals are proved for a flux function of the 
form H(x,t,p) = f(p) + F(x,t) under the periodic setting, with many details of large 
time behaviors of the Lax-Friedrichs scheme. Still it seems very hard to obtain results 
similar to [13] for more general flux functions by the standard approach based on the 
L 1 -framework. 

Recently the author announced a stochastic and variational approach to the Lax- 
Friedrichs scheme [18], where stability and convergence are proved based on the low of 
large numbers in hyperbolic scaling limit of random walks and the calculus of variations 
in the theory of viscosity solutions of the Hamilton- Jacobi equations with constants c 
and h(c) 

(1.2) v t + H(x,t,c + v x ) = h(c). 

This is a finite difference version of the work [9] on the stochastic and variational approach 
to the vanishing viscosity method. Now we briefly refer to our stochastic and variational 
approach. Consider initial value problems of the inviscid hyperbolic scalar conservation 
law and the corresponding Hamilton- Jacobi equation 



1.3) 



u t + H(x,t,c + u) x = in Tx (0,7*|, 

u(x,0) = u°(x) G L°°(T) on T, [u°(x)dx = 0, || w° |U-< r, 

Jr 

v t + H{x,t,c + v x ) = h(c) in Tx (0,!T], 
v(x, 0) = v°(x) G Lip(T) on T, || v° ||l°°< r, 



(Al) H(x,t,p) : T 2 x R^R, C 2 (A2) H pp > (A3) lim ,1, = +oo. 



where c is a parameter varying within an interval [cq, C\], T := R/Z is the standard 
torus, h(c) is a continuous function and r > is a constant. We arbitrarily fix T, r and 
[c ,ci]. Note that (11.31) and (11.41) is equivalent in the sense that the entropy solution u 
or viscosity solution v is derived from the other if u° = v®. In particular we have u = v x 
(see e.g. [I]). We always assume that u° = v®. Our flux function H is assumed to satisfy 
the following (Al)-(A4): 

H(x,t ,p) 

|pH+oo \p\ 

By (Al)-(A3), we have the Legendre transform L(x, t, £) of H(x, t, •), which is now given 
by 

L(x,t,£) = sup{£p- H(x,t,p)} 

PGR 

and satisfies 

(Al)' L(x, t,():T 2 xl^M, C 2 (A2)' L« > (A3)' lim ^ll^l = +00 . 

ICK+oo |^| 

The last assumption is 

(A4) There exists a > such that \L X \ < oc(\L\ + 1). 
We discretize the equation of ( 11. 3D by the Lax-Friedrichs scheme as follows: 



k+l ("m+"m + 2) 



H(x m+2 , t k , c + u k m+2 ) - H(x m , t k , c + u^ 



'1 5) m +l 2 _|_ " \~m-r^i -K-i " ' ~"m+'i/ ~~ \~nii -ki - i ~~ra) q 
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We also discretize the equation of ( II .4p by the following scheme 



v 



fc+i 



m+l 
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(1.6) 
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+ H(x m ,t k ,c + 



m—1 



) = Kc). 



At 



Note that ( 11.51) and fl 1.6ft are also equivalent in the sense that or v^ +1 is derived from 
the other. In particular we have 



which is an important relation in this paper. In the stochastic and variational approach, 
"stochastic" comes from the numerical viscosity intrinsic to (II. 5p and (II. 6p . and "vari- 
ational" from the variational structures of Hamilton- Jacobi equations. We state several 
points of the stochastic and variational approach [T8] : 

(1) Stochastic and variational representation formulas (value functions) for v^+i and 
v^ m are given. 

(2) Stability of the Lax-Friedrichs scheme up to arbitrary T > is derived from vari- 
ational techniques. 

(3) Almost everywhere pointwise convergence of to u = v x is proved. In particular 
this yields the uniform convergence except neighborhoods of shocks with arbitrarily 
small measure. 



(4) Uniform convergence of v^ +1 to v with an error O(VAx) is proved from a stochastic 




(5) Random walks play a role as "characteristic curves" of the difference equations, 
which converge to the genuine ones. 

The purpose of this paper is to show further results on the Lax-Friedrichs scheme, 
under the assumptions (Al)-(A4), based on (l)-(5) and with techniques of the theory of 
viscosity solutions of Hamilton- Jacobi equations. Although our arguments here are on 
the Lax-Friedrichs scheme, other finite difference schemes with the numerical viscosity 
are available as well. The main results are on 

1. Time global stability of the Lax-Friedrichs scheme with a fixed meshed size. 

2. Error estimates for entropy solutions. 

It is proved that genuine entropy solutions at t — 1 are uniformly bounded, which 
is independent of the magnitude of initial datas. Since difference entropy solutions are 
"close" to the genuine ones for small mesh size, difference entropy solutions at t — 1 
are also uniformly bounded. Due to the periodic setting, iteration of time-1 analysis 
yields time global properties. Combining these facts, we obtain time global stability 
of the Lax-Friedrichs scheme. As a result, we can show the long time behavior of the 
Lax-Friedrichs scheme that any solutions associated with c fall into the time periodic 
state uniquely determined by each c, which means that we obtain the unique space- 
time periodic difference entropy solution and the unique, up to constant, space-time 
periodic difference viscosity solution for each c. They approximate a genuine Z 2 -periodic 




m—1 




viewpoint. 



3 



entropy/ viscosity solution of fll.ip /( !L2|) . In the periodic states, we naturally have the 
notion of the effective Hamiltonian for difference Hamilton- Jacobi equation ( II. 6p . We 
also prove convergence with an error estimates of the effective Hamiltonian to the genuine 
one, showing its properties. 

It is known that the optimal L 1 -error estimate between u k m and u is 0{\fKx) in 
the case of H(x,t,p) = H(p) [15] . The upper bound 0{\fKx) is due to properties of 
functions of bounded variation [TTJ. It is not clear whether [TTJ is applicable to the case 
of our general flux functions. We show, with different approach, an L 1 -error estimate of 
O(Axi). Our error estimate is based on 0(V Ax) which arises as an "error" between 
random walks for the Lax-Friedrichs scheme and their space-time continuous limit under 
hyperbolic scaling (that is a backward characteristic curve). For a technical reason, we 
further lose 1/4 in the case of general flux function H(x, t,p). In addition we show that, 
if the genuine entropy solution is Lipschitz, a C°-error estimate of O(Axi) is available. 

Unlike the case of initial value problems, it is challenging to show convergence of full 
sequences and an error estimate for Z 2 -periodic entropy/ viscosity solutions of (11.11) / (II. 2p . 
because uniqueness of genuine Z 2 -periodic solutions with respect to c is not valid in 
general. However we can manage the special case where a genuine Z 2 -periodic entropy 
solution u with some c is C l and the dynamics of its characteristic curves C*(s) := (q(s) 
mod 1, s mod 1) are C 1 -conjugate to that of the linear flow on T 2 with a Diophantine 
rotation vector. Existence of such a solution u is known as a KAM torus in Hamiltonian 
dynamics (see e.g. [12], [16], pjj]). We show a C°-error estimate depending on the 
Diophantine condition of the rotation vector, which is a result on approximation of 
KAM tori. Our proof is based on the fact that one orbit of the linear flow on T 2 with a 
Diophantine rotation vector is ergodic on T 2 . It may be impossible to obtain the same 
result without our stochastic and variational approach to the Lax-Friedrichs scheme, 
because if so there is no way to track the genuine characteristic curves, namely the 
result relies on the fact that random walks are close to the genuine ergodic characteristic 
curves. 

Finally we state that our motivation comes from not only continuum mechanics but 
also a recent theory of Lagrangian and Hamiltonian dynamics called the Aubry-Mather 
theory or the weak KAM theory (also from the classical KAM theory) [SJ, [HI], [7]. Our 
periodic setting is standard and Z 2 -periodic entropy/ viscosity solutions of (II. ip / (11.21) and 
the effective Hamiltonian play central roles in the weak KAM theory. The results of this 
paper provide basic tools of numerical analysis of the weak KAM theory through finite 
difference approximation. We remark that it is better to approximate entropy solutions 
and characteristic curves as well as viscosity solutions from the view point of accuracy, 
because the central objects in the weak KAM theory such as KAM tori, Aubry-Mather 
sets, effective Hamiltonians, calibrated curves etc. are derived from the derivative of 
viscosity solutions or entropy solutions. The "derivative" of numerical viscosity solutions 
obtained through a scheme which has no relation to entropy solutions is not accurate 
in general. Some efforts in finite difference approximation to the weak KAM theory are 
found in [13]. However due to the absence of the stochastic and variational approach 
to the Lax-Friedrichs scheme, the results are mathematically restricted. We also point 
out [10], [2] for results on smooth approximation methods of the weak KAM theory 
by the vanishing viscosity method. In particular [2] successfully exploits the stochastic 
and variational approach to the vanishing viscosity method given in [5] , where stochastic 
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ODEs with the standard Brownian motion approximate the genuine characteristic curves. 
Our results allow ones to develop more on finite difference approximation methods of 
the weak KAM theory. 

2 Preliminary Result 

We state several important background or preliminary results. 
2.1 Entropy Solution and Viscosity Solution 

It is well-known that the viscosity solution v of fll.4p . which is Lipschitz, is characterized 
by calculus of variations: The value of v at each point (x,t) G T x (0, T], T £ (0, oo) is 
given by 

(2.1) v(x,t) = inf ( rL c (7( S ), S ,7 / (^))^ + ^°(7(0))) + ^(cK 

where AC is the family of absolutely continuous curves 7 : [0, t] — > T and 

L c (x,t,0 :=L(xM)-c£ 

is the Legendre transform of H(x,t,c+ •). We can find a minimizing curve 7* of (12. ip . 
which is a backward characteristic curve of (II. ip and (11.21) . and a C 2 -solution of the 
Euler-Lagrange equation generated by the Lagrangian L c . On each minimizing curve, v 
is differentiable with respect to x: 

(2.2) v x (j*(s),s) = Ll( 1 *(s),s, 1 *'(s)) for < s < t. 

We say that a point (x,t) is a regular point of v or just regular, if there exists v x (x, t). 
Since v is Lipschitz, almost every points are regular. In particular, if (x, t) is regular, a 
minimizing curve 7* for (12.11) is unique and (12.21) holds for s = t. 

Usually the entropy solution u of (11.31) is defined as an element of C°((0, T]; L X (T)). 
Here we always take the representative element of such u given by v x , which is still 
denoted by u. If (x,t) is regular and 7* is the unique minimizing curve for v(x,t), the 
value of the entropy solution u = v x at the point (x, t) is given by 

u(x, t) = f L c x ( 7 *(s), s, 7 *'(s))ds + u (7*(0)), 
J 

where u° is supposed to be rarefaction- free, or equivalently v is semiconcave. Otherwise 
u°(7*(0)) needs to be replaced with L|(7*(0), 0, 7*'(0)). In particular we have for any 

T£[0,t) 

u(x, t) = J' Ll(>f(s), s, r(s))ds + ^( 7 *(r), r, 7 *'(r)). 

For more details, see e.g. [I], [1]. 

We introduce the solution operators of (11.31) and (11.41) 

0* : L^ (T) 3 vP 1 — y u(-,t) £ L°°(T), ^ : Lip r (T) 3 v° h-> v(-,t) £ Lip(T), 

5 



where L£° (T) is the set of all functions u° G L°°(T) with || u° r and J T u°dx = 0, 

and Lip r (T) is the set of all Lipschitz functions on T with a Lipschitz constant bounded 
by r. When we specify the value of c, we write </>*(•; c), ^(-; c), u c , v c . 

We would like to prove an a priori boundedness of u(x, t) = v x (x, t). 

Proposition 2.1. For each t G (0, T], there exists a constant fii{t) > independent of 
r, c G [co,Ci] and initial datas v ,u° for which we have 

II 0V; c) || LOO < || c), || L ~< /?i(t). 

Proof. Fix t G (0, T]. Due to (12. 2p . which holds for s = t if (x, t) is regular, it is enough 
to estimate L^(7*(t), t, 7*'(£)) for each minimizing curve 7*. We prepare two lemmas. 

Lemma 2.2. Let 7* &e a minimizing curve for v(x,t). Set y := 7*(0). Then 7* attains 

-t 



inf / L c ( 7 ( S ), S ,Y(s))^. 



7eAC,7(t)=a;,7(0)=?/ J 

Proof. If not, there exists 7" such that 

L C ( 7 «( S ), s, 7 fl '(^))^ < f L c (j*(s), s, j*'(s))ds. 

Jo 

Since ^°(7 t) (0)) = y = w°(7*(0)), we have 

L c (^(s),s,^'(s))ds + v°(^(0)) < [ L c ( 1 *{s),s,'y*'{s))ds + v { 1 *(0)). 

Jo 

Therefore 7* is not a minimizing curve for v(x,t), which is a contradiction. □ 
Define the set 

T(t) := (7° I 7 C attains inf [ L c (j(s), s, j'(s))ds } 2, y G T, c G [c , Ci] 1 . 

By Lemma 12.21 any minimizing curve 7* for v(x,t), x G T belongs to (To be 

precise, we should take 7* mod 1, but this is not important due to the periodic setting.) 

Lemma 2.3. 1. There exists a constant C\(t) > such that for any x, y G T we have 
a C 1 - curve 7 which satisfies 

y{t)=x, 7(0) =2/, f L c {i{s),s,i{s))ds<Ci(t). 

Jo 

In Particular, any 7 C G T(t) satisfies 

L c (l c (s),s,^f c/ (s))ds < Ci(t) 



2. There exists a constant C^it) > such that for any 7 C G T(t) we have r G (0,t) 
which satisfies 

h c '(r)\<C 2 (t). 
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3. There exists a constant C$(t) > such that for any 7 C G T(t) we have 

|^( 7 C ( S ), S , 7 C '( S ))| <C 3 (t), se[0,t}. 

Proof. 1. Consider 7(5) := x + — t). Since \x — y\ < 1, we have | 7 '(s)| < t~ x . 

Therefore we obtain 

/ L c ( 7 ( S ), S ,7'( S ))ds < sup |L c (x,s,0|t 

JO x,s£T,\^\<t- 1 ,ce[c ,ci] 

Set Ci(t) : = sup a , iSeT! |^|< t -i )C6 [ COjCl ] \L c (x, s, £)\t and 1. is proved. 

2. We have r G (0, t) which satisfies, due to 1. and the minimizing property of 7 C , 

Cx(t) > / L C ( 7 C ( S ), a, 7 c '(s))^ = L c ( 7 c (r), r, 7 c '(r))t. 
•/ 

By (A3), | 7 c/ ( r )| must be bounded by a constant C^(£) independent of 7 C G T(t). 

3. Note that 7 C is a C 2 -solution of the Euler-Lagrange equation generated by L c 

|^(7 c ( S ), S , 7 c/ ( 5 )) = ^(7 c ( S ), S ,7 c/ (^))- 

It follows from (Al)-(A4) that there exists a\ for which we have \L C X \ < ai(\L c \ + 1) 
for any c G [c ,Ci], and that := | min{0, inf L c }\ is bounded. We have r* G [0, t] 

which attains the maximum of |L^( 7 c (s), s,7 c/ (s))|. Suppose that r* ^ r, where r is the 
number in 2. 

^(7 c (sU 7 c/ (s))rf S | = |^(7 c (r*),r*,7 c/ (r*))-^( 7 c (r),r,7 c/ (r))| 



< r iw),^w)i<fa 
«/ 

< r ai (l + |L c (7 c (s), S ,7 c/ ( S ))|)d S 
>/ 

< ai / 1 + (L C ( 7 C ( S ), S ,7 C/ ( S )) + L,) + L^s 

Jo 



= ai(2L + + l)t + «! / L c (7 c (s), S , 7 c/ ( S ))ds 
< a 1 (2L* + l)t + a 1 C 1 (t). 



Therefore, setting 



C 3 (t) :=ai(2L» + l)t + aiCi(t)+ sup |L^(x,s,OI, 

x,seT,|5|<c 2 (t),ce[co, Cl ] 

we obtain for < s < t 

mms)^,^'^))] < |^(7 c (t*),t*,7 c/ (t*))| <C 3 (t). 
The case where r* = r is included by the above inequality. □ 
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Since (x,t) is regular for almost every x G T with each fixed t and v x (x,t) = 
L|(7*(t), t, J*'(t)) holds almost every x G T, we obtain Proposition 12.11 by setting 
Pi{t):=C 3 (t). □ 

We show continuity of </>*(f°; c) and ^(vq] c) with respect to u° and c. 

Proposition 2.4. Fzx £ G (0,T]. For each sequence v® —> v° uniformly and c J — > c as 
j — > oo (v® x is not necessarily convergent), we have 

^(Vj] c 7 ) -» ip^v ; c) uniformly, 4>\v jx ] c 7 ) -»■ c) in L X (T) as j ->■ oo. 

Proof. By the variational representation, we have 

^(v°;c)(x) = f L( 1 *(s),s, 1 * l (s))-c 1 *'(s)ds + v°( 1 *(0)) + h(c)t, 
Jo 

^tyVX*) = [ L(^(s),s,j*'(s))-cJ 1 *'(s)ds + w?(tJ(0)) + h(C> 





-c)Y'(s)ds + v° j (Y(0)) 


-«V(o)) 




- h(c))t, 




Jo 


-c) 1 ;'(s)ds + v^ 1 *(0)) 


-«°(^(0)) 




-h{c))t. 





and hence 

^(ujjc 7 ')^) -^(uVXz) < 

V>*(uJ; cOfa) -^*(u°;c)(z) > 



It follows from 3. of Lemma, 12.31 that any minimizing curves for v(x,t) are Lipschitz 
with a common Lipschitz constant for all x G T and v° G Lip r (T). Since h is continuous, 
we conclude ^(v®; e 7 ) — >■ ^(v ; c) uniformly as j — > oo. 

Let x G T be a common regular point of all ^(v®; e 7 ), j = 1, 2, 3, • • • . Almost every 
point is such. Due to a variational technique, we see that 7* — >■ 7* uniformly and 
7*' — > 7*' in L? as j 00 (see e.g. [18], Lemma 3.4). Note that for each < r < t 

ft 



<P\vl- c){x) = iP\v»; c) x (x) = J L x (j*(s), s, 7>))<fe + W(r), r, T *'(r)) - c, 
<P{*U c J )(x) = <t) x {x) = f L„(tJ(«), «, < (a))tfa + ^(tJ(t), r, t$'(t)) - c 7 . 



For any £ > 0, there exists J such that, if j > J, we have || 7* — 7* ||c°< £ > II 7j' — 7*' |U 2 < 
ey^- Note that we have r depending on j > J such that |7*'( r ) — 7*'( r )l — e - Therefore 
we conclude 0*(f° x ; e 7 ) — >■ °; c) pointwise almost everywhere. This immediately leads 
to L 1 (T)-convergence. □ 



2.2 Stochastic and Variational Approach to the Lax-Friedrichs 
scheme 

We state several results of the stochastic and variational approach to the Lax-Friedrichs 
scheme shown in [18]. Let N, K be natural numbers with N < K. The mesh size 
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A = (Ax, At) is defined by Ax := (2N)^ and At := (2K)~ l . Set A := At/Ax, 
x m := mAx formal and t k := kAt for k = 0, 1, 2, • • • . For x G M. and £ > 0, the 
notation m(x), k(t) denote the integers m, k for which x G [x m , x m +2Ax),t G \t k , t^+At). 
Let (Axli) x (AtZ> ) be the set of all (x m ,t k ) and 

Q even C (AxZ) x (AtZ> ), Qodd C (AxZ) x (AtZ> ) 

be the set of all (x m , t k ) with k = 0, 1, 2, • • • and m G Z such that m + k =even, odd. 
We call Q even , Qodd the even-grid, odd-grid. We consider the discretization of (jl.3p by 
the Lax-Freidrichs scheme in Q even : 

U m+l 2 ~ H(x m j r 2, tk, C + M^ +2 ) ~~ H(x m , tf~, C + M^) 



(2.3) ^ 



At 2Ax 



m±2N ~ u mi 



where for m = even 

fX m +Ax 



(2.4) u° A (x) := — — / u°(y)dy for x G [x m — Ax, x m + Ax). 

ZAX J Xrn -Ax 

Note that • 2Ax is conservative with respect to k and is zero for 

{m | 0<m<2N, m+k=even\ 

u° which has the zero-average. We also discretize ( 11.41) in Q m'- 



(2.5) 



v k+l _ ("m-l+^m+l) k _ k 

^ + ( W fe , c + m+1 2A ^- 1 ) = h(c) 



V m+1 — v A( X rn+l), V m+1±2N ~ V : 



m+li 



where, in addition to the assumption u° = v®, we assume that 

PX 

(2.6)v° A (x) := v°(-Ax) + / u° A (y)dy (v° A (x m+1 ) = v°(x m+1 ) for m =even). 



Ax 



Note that u° A — > u° in L X (T) and v A — > u° uniformly with || v A — v° u ° -2Ax, 

as A 0. We introduce the difference operators: 

fc+i _ <-i+»4+i k _ k 

2Wj+i : =I^ — -J? , D x w k m+1 := Wm+1 ™ m ~\ 

At m+1 2Ax 

The two problems (12.31) and (12.51) are equivalent under (12.41) and (12.61) . In particular we 
have -D x i^j +1 = u k m p2]. Let «a be the step function derived from the solution u k m of 
(12.31) . namely 

u A (x,t) := u* for (x,t) G [x m _i,x m+ i) x [t k ,t k+1 ). 

Let f a be the linear interpolation with respect to the space variable derived from the 
solution v^ +1 of (I2.5p . namely 

v A (x, t) := v k m _ x + D x v k n+1 ■ (x - x m _i) for (x, t) G [x m -i, x m +i) * [tk, t k +i)- 
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We remark that v&(x, •) is a step function for each fixed x and that (v&) x = ma- 

We introduce space-time inhomogeneous random walks in Q m, which correspond to 
characteristic curves of (II. 3p and (jl.4p . For each point (x n ,t [+1 ) G Qodd, we introduce 
backward random walks 7 which start from x n at and move by ± Ax in each backward 
time step: 

7 = { 7 fc } fc=0iW+1 , 7 l+1 = x n , 1 k+1 - 1 k = ±Ax. 
More precisely, we introduce the following for each (x n ,ti + i) G Qodd- 

X k := {x m I (x m ,t k ) G Q dd, \x m -x n \ < (I + 1 - A;) Ax} for fc < I + 1, 
G := |J x {t fe }) C Q dd, 

\<k<l+\ 

i:G3 (x m ,t k ) ^ g [-A -1 , A" 1 ], A = At/Ax, 

11 

p : G 3 (x m ,t k ) H- p^ := - - -A^ G [0, 1], 
P ■ G 3 (x m , t fc ) 1 ^ p^ := - + -A£* G [0, 1], 

7 : {0, 1, 2, • • • , / + 1} 3 k H> 7 fc G 7 i+1 = x n , 7 fc+1 - 7 fc = ±Ax, 
f2 : family of the above 7. 

We regard p^ (respectively pjjj as a transition probability from (x m , to (x m + Ax, t k — 
At) (from (x m ,tk) to (x m — Ai,4 — At)). Note that £ is a control of random walks, 
which plays a role of a velocity field on the grid. We define the density of each path 
7 G Q as 

M7):= II P(7 fc ; 7 fc " 1 ), 
i<fc<«+i 

where p(7 fc ,7 fc_1 ) = P^( 7 fc) (respectively p^^/tO if 7 fc — 7 fc_1 = —Ax (Ax). The density 
//(•) = //(•;£) yields a probability measure of f2, namely 

prob(A) = ^2 ^(7; f° r Acfi. 

The expectation with respect to this probability measure is denoted by E^.-^, namely 
for a random variable / : Q — > K. 

^(•;0[/(7)] :=^p:(7;0/(7)- 

We use the notation 7 as the symbol of random walks or a sample path. If necessary, 
we write 7 = 7(x n , t [+1 ; £) in order to specify its initial point and control. 

We state an important result on scaling limit of inhomogeneous random walks. Let 
vil) = {v k (l)}k=o,i,2,-,i+i, 7 G f2 be a random variable which is induced by a random 
walk 7 = 7(x n , t[ + i, £) and is defined as 

V l+1 : =1 l+ \ r / fc ( 7 ):= 7 m - £ M At for < k < 1 

k<k'<l+l 
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Proposition 2.5. ({Hj) Set d k := E^,^ - ^( 7 )| 2 ] and d k := E^,^ - ^( 7 )|] 
for < k < I + 1 . T/ien we have 

4.1+1 _ J-k 

(d k f <a k < Ax. 

A 

If we take hyperbolic scaling limit, namely A = (Ax, At) — > under 

< A < A = At/ Ax < Ax, 

d k and yd* always tend to with the order 0(a/Ax). Note that the variance does not 
necessarily do so for inhomogeneous random walks. We refer to p2] for more details 
of hyperbolic scaling limit of inhomogeneous random walks. We always take the limit 
A — > under hyperbolic scaling. 

Now we state results on the stochastic and variational approach to the Lax-Friedrichs 
scheme: 

Theorem 2.6 ([IB]). There exists Ai > (depending on T, [c , c±] and r) such that for 
any small A = (Ax, At) with A = At/ Ax < Ai we have the following: 

1. The following expectation given by 7 = j(x n , tt+i] £) 

s M , e) [ ^ c (7 fc ,4-i,eV)) At+v A(7 )] +h( C )t l+1 

0<fc</+l 

has the infimum with respect to £ : G — >■ [— A _1 ,A _1 ] /or eac/i /ixerf n G Z and 
< / + 1 < fc(T). T7ie infimum is attained by £* £/iat satisfies |£*| < A7 1 < A -1 . 

The solution v 1 ^ 1 of A2.5\) satisfies for each nfZ and < / + 1 < k(T) 

v l n +1 = mi E K , [ £ ^ c (7 fc ,4-i,e (7fc) )At + ^( 7 )] +h(c)t l+l . 

0<fc<Z+l 

3. The minimizing velocity field £* for each v 1 ^ 1 is unique and satisfies in G 

E^(x m ,tf.,^ m ) = D x V m+l (O- £ m = H p (x m , tfc, c + D x v m+1 )). 

4- LetS,* be the minimizing velocity field for v 1 ^ 1 ■ Let 7 = 7 (x n , ti+i, £*) andfi(-;C,*) be 
the minimizing random walk and its probability measure. Let £* be the minimizing 
velocity field for t^+V Let 7 = 7 (x n+ 2, tz+i; £*) and /i(-;£*) fre £/ie minimizing 
random walk and its probability measure. Then u l ^l\ = D x v l r ^ 2 satisfies 

u l + + \ < E K , e) [ L:( 1 k ,t k . 1 ,C k m{ , k) )At + u A ( 1 + Ax)]+O(Ax), 

0<k<l+l 

> £ ^(7 fc ,Vi,lt (W )At + <(7°-Ax)]+0(Ax), 

0<fe</+l 

where O(Ax) stands for a number of (—9 Ax, 9 Ax) with 9 > independent of Ax. 
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Now we take hyperbolic scaling limit: 

5. Let v be the viscosity solution of jl-4\) - Then for each t G [0,T] 

v A{',t) — > v(-,t) uniformly on T as A — >■ 0. 

In particular, we have an error estimate: There exists fii > independent of A, 
c G [cq, Ci] and initial datas v° G Lip r (T) such that 

SUp || V A (-,t) - v(-,t) ||c"(T)< fcvAx- 
t£[0,T] 

6. Let (x, t) G T x (0, T] be a regular point and 7* : [0, t] — > R be the minimizing curve 
for v(x, t). Let (x n , be a point of [x — 2 Ax, x + 2 Ax) X [t — At, t + At) and 
7a : [0,£] — )■ M 6e t/ie linear interpolation of the random walk 7 = 7(x n , £*) 
given 6?/ t/ie minimizing velocity field /or . Tnen 

7a — >■ 7* uniformly on [0,t] in probability as A — > 0. 

In particular, the average 0/7A converges uniformly to 7* as A 0. 

7. Let u = v x be the entropy solution of U.3\) . Then for each regular point (x,t) G 
T x [0, T] 

u A (x,t) — )■ u(x,t) as A 0. 

In particular, ma converges uniformly to u on (T x [0, T]) \ 0, where is a neigh- 
borhood of the set of points of discontinuity of u with an arbitrarily small measure. 

Note that 1. and 3. give the stability condition of the Lax-Friedrichs scheme, which 
is called the CFL- condition 

\\H p (x m ,t k ,c + u k m ))\ < 1. 

We further state preliminary results on the Lax-Friedrichs scheme. Introduce the 
solution operators of (12.31) and (I2.5P 

</> A : L^ (T) 3 u° ^ u A (-,t) G L°°(T), Va : Li Pr (T) 3 v° ^ v A {-,t) G Lip{T). 

When we specify the value of c, we write 4> A (-] c), V'a(') c )> u % u m( c )i v< ki v m+i( c )- Note 
that we first get the step function u° A from u° as (12.41) and then maps u° A to u A (-,t) by 
A . Similarly we first get the piecewise linear function -o A from v° as (12. 6p which satisfies 
u° = v® and then maps v A to v A (-, t) by i/) A . 

Proposition 2.7. Fix t G [0, T]. For each sequence v*j — > u° uniformly and cP — )■ c as 

j — J- 00 zs not necessarily convergent), we have 

■0 A (v°; c 3 ') ->■ ^a( w °; c ) uniformly, <p A {v° jx ; c?) -» A (f°; c) in L l (Y) as j -> 00. 
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Proof. It is enough to show c?)(x n ) —> c)(x n ) uniformly with respect 

to x n as j — >■ oo. Using the stochastic and variational representation, we have 

^(v°;c)(x n ) = £ M , e) [ ^(7 fc ,Vi,et (7fe) )-cr^)At + <( T )" 

0<fc<Z+l 

+/i(c)t z+ i, 

V4 +1 KV')(^) = £ ^(7 fc ,4-i,e j t {7 ,))-c J e^ (7fc )At + ^ A (7 

0<fc<i+l 

+h(d>)t l+1 , 



where are minimizing velocity field. Hence, by the stochastic and variational 

representation again, we have 

VlTfyVXsn) - ^l +1 (v°;c)(x n ) 

< E -(^-^t (7 ,)At + ^ A ( 7 °)-<( 7 ) 

0<fc<«+l 

+(/i(c J ")-Mc))tj + i, 
^ +1 (vJ;c*)M - ^ A +1 (^°;c)(x n ) 

> E ^q)[ E -(^-^ 7 *) A * + w iA(7°)-«i(7°) 

0<fc<«+l 

+(/iry)-/i( c ))t I+1 . 

Since are uniformly bounded, we conclude the assertion. 

The second convergence follows from the first one and the following relation: 

, t f o n / \ 0? ; c J ' ) (x m+ i ) - V> A (vj ; c j ) (x m _i ) 
4>\{.v u jx ]c , ){x m ) = J - 



2Ax 



□ 



We show details of the one-sided Lipschitz condition of it^ which leads to the entropy 
condition of u. We remake that the one-sided Lipschitz condition of u k m is not necessary 
for Theorem \2.6] namely u k m converges to the exact entropy solution regardless of the 
condition. However we want to investigate the one-sided Lipschitz condition in addition 
to Theorem 12.61 in order to obtain further results. It is easy to prove the one-sided 
Lipschitz condition through v^ +1 to be a "semiconcave function" of a finite difference 
version, using the stochastic and variational representation. Unlike the case of exact 
viscosity solutions, however, we need the assumption of v to be semiconcave. Even if 
v° is semiconcave, it is not easy to estimate time-dependency of the one-sided Lipschitz 
constant. Therefore we use a direct method similar to that of Lemma 2 in [T3]. The 
direct method is available up to an arbitrary T > 0, because we already know about 
boundedness of difference solutions up to T due to Theorem 12. 61 Introduce the following 
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notation with Ai in Theorem 12.61 

u* := sup \L^(x,t,^)\ (note that |it*J <u*), 

x.teT^elco.cil.l^l^Aj" 1 

H* x := sup c + m)|, H* xp : = sup |if xp (x, t, c + 

x,tST,c£[co,ci] ,|u|<u* a:,tST,cG[co,ci],|u|<ii* 

#* := inf \Hpp(x,t, c + u)\ (H* > due to (A2)), 

x,teJ,ce[co,ci],\u\<u* 



1 2H* / /if* \ 2 Off 

V ■= max{2^ p + -H* pp + i^}, E * := ^ + M-g) + — 



-pp V pp pp 



Proposition 2.8. Let Ai > be that of Theorem \2.6l Suppose that A = (Ax, At) 

I* ) 

ppj 



satisfies A = At/ Ax < A 1; At < min{(2r/)- 1 , (2E*H*)~ 1 } ; 



(2.7) sup X(\H p (x,t,c + u)\ + H* xp -2Ax) < 1, A < — -. 

s,teT,c6[co,ci],|u|<u* r pp 

Then the following holds: 

1. We have for 1 < k < k{T) 

E k A := max ^ +2 ~^ <^i (t k = kAt). 
m 2Ax ~ H; p t k ' 7 

2. < £7*, we /iawe < E* for 1 < k < k{T). 

2e 2 n 

3. Ifk> k(r]- 1 ), we have E k A < —f. 

4- If ' u k m is extended to k — > oo with < w* ; we have limsupi?^ < E*. 

k~ >oo 

Proof. Set 2;^ := — M m- Using the difference equation and Taylor's formula, we 
have 

y k i „fc 



2 m+l = o 2~Ar ^^ m+4 ' ^' ° Mm+4 '' _ ^( X m+2, tk, C + 1t m _ 



f H (x m+2 , t fc , c + u k m+i ) - H(x m+2 , t k , c + u k m+2 ) 
f H(x m , t k , c + u k m ) - H(x m+2 , t k , c + u k m ) 
^H(x m+2 , t k , c + u k m )- H(x m+2 , t k , c + u k m+2 )} 

,~ + —H p (x rn +2, t k , C + M^))^ + ~~ -^Hp{ x rn+2, tk, C + !iJ, +2 ))4 +2 



At 



2Ax 



{(iJ x .(x m+ 2 5 t k , C + M m+4 ) - H x (x m+2 ,t k ,c + u k n ))(2Ax) 

+ \h pp ■ {z k m+ 2? + \H„ ■ (z k m f + \H XX ■ {2Axf + l -H xx ■ {2Axf} 



{\ + ^H P (x m +2, t k , c + u k m+2 ) - ^H xp ■ 2Ax)z k m 
+{ 2 — ~^H p (x m+ 2, t k , c + u k mJv2 ) — —H xp ■ 2Ax}z^ n+2 
^{\ H pp ■ (^ +2 ) 2 + \h„ ■ (z k J 2 + \h xx ■ {2Axf + \h xx • (2Ax) 2 } 
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By the first inequality in (12. 7p . it holds that 

{\ ± ^H p (x m+2 , t k , c + u k m+2 ) - ^H xp ■ 2Ax} > 0. 
Hence, setting z^ := max{^,^}, we obtain 



H* At 

z k m \\ < {l-2H xp At)^ + Hl x .2AxAt-^ — {i k m )\ 



~k 17* ~k 

*m+l ^ otj \,\ ^rn , tt* a + PP - ■' ^ 



< (l-2H xp At)-f- + H: x At-^At(, 



2Ax ~ v p '2Ax xx 2 y 2Ax' 

Note that g(y) := (1 - 2H xp At)y + H xx At - (^-At)y 2 is monotone increasing, if y < 
(1 - 2H xp At)/(H* p At). It follows from the second inequality in ( TP) and At < (2r/)~ 1 
that E\ < r/2Ax < (1 -2H* xp At) / (H* pp At) < (1 - 2H xp At) / (H* p At) for all initial datas. 
Suppose that E k A < (1 - 2H* xp At) /(H* pp At) . Then we obtain 

(2.8) E k A +l <E k A - At{^(E k A ) 2 - 2H* xp E k A - H* xx }. 

Since At < (2E*H; p )~\ we have E* < (1 - 2H* xp At)/(H; p At). If < £| < (1 - 
2H* p At)/(H; p At), we have < E k A . li E k A < E* , we have > E\ but < 
Therefore it holds that E k A < (1 - 2H* xp At) / {H* pp At) for all < fc < fc(T). Thus (l2~8l) 
holds for all < k < k(T). We can check that E A may increase up to E*, if E A < E*, 
and i?^ are bounded by a monotone decreasing sequence with lim sup^^ E A < E*, if 
£^ > E*. 

Now we follow Lemma 2 in [T3]. Set V k := i?^ + 1 > 1. Then we have 

H* 

V k+1 < (1 + 7]At)V k - ~^{V k ) 2 
Set W k := (1 - rjAt) k V k , k > 0. Then we have for fc > 1 

< (1 - V At)(l + V At)W k - -J^At(W k ) 2 (l -r]At)~ k+1 

< W k _ ?h.At(W k ) 2 . 

Consider w'(t) = -^f (w(t)) 2 , w(0) = w° := 2/(H* pp At). The solution satisfies 



= if T - 

PP J- i 1 



2 



H* t 



We can show that W k < w(kAt) for k > 1 as follows: Note that w(At) = l/(H* p At) 
and W/ 1 = (1 - r]At)(E A + 1) < r/(2Ax) + 1. It follows from At < (2T/)" 1 and the 
second inequality in (12. 7\i that W 1 < w(At). Suppose that W k < w(kAt) for some 

H * A t: 

k > 1. Then, since g(y) := y ^ — y 2 is monotone increasing for y < At), 

io(fcAt) < l/(H* p At) and u/' > 0, we have 

= iw(A;At) + Atw'(kAt) = w(kAt + At) - -w" ■ (At) 2 



2 



< w({k + l)At). 
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Thus we obtain 

„h a x k 2 / . v vkAt 2 2e r * th 1 

E ^ (1 " " At) 5^7 £ f 1 " " Atr w 7^ £ ^ 

Set f(t) := jfrj- The minimum of / is /(i]^ 1 ) = IS 2 -, which is larger than Since 
E\ with £ A > are bounded by monotone decreasing numbers, we conclude that 
Ei < f{ri~ l + At) < §^ for k > k(^ 1 ). □ 



3 Time Global Stability and Long Time Behavior 

We prove time global stability of the Lax-Friedrichs scheme with a fixed mesh size. Then 
we show long time behaviors of the scheme in which each difference solution fall into a 
time periodic state with the period 1, obtaining space-time periodic difference solutions. 
The notion of the effective Hamiltonian of (11.61) arises. 



3.1 Time Global Stability 

The main result of this section is the following: 

Theorem 3.1. There exists \\ > and 5 > such that, if A = (Ax, At) satisfies 
< Ao < A := At/ Ax < Ai and Ax < 5, the Lax-Friedrichs scheme starting from any 
u° G L£° (T) works up to an arbitrary time index with the CFL-condition 

\H p (x m , tk, c + u^JI < A^ 1 < A -1 for all m G Z and k G Z + . 

In order to prove this theorem, we need uniform boundedness of || </> A (u°;c) ||x,oo with 
respect to (u°; c) similar to Proposition 12.11 First we observe the following lemma: 

Lemma 3.2. Let Ai > be that of Theorem \2.6\ and A = (Ax, At) be such that < 
Ao < A = At/ Ax < Ai. Fix arbitrarily t G (0,T]. Then for any e > there exists 
5 = 5(e; t) > such that, if Ax < 5, we have 

sup || A (w°; c) - c) || L i(T)< e. 

« 6L^ (T),c6[c ,ci] 

Proof. If not, for some Eq > and Sj — > 0, we have Axj < 5j such that 
(3.1) sup || 4>V{u°; c) - 4>\u°; c) \\ L ^r)> e , 

w°6i~ (T),ce[co,ci] 

where Aj = (Axj, XAxj). We show that, for each j, there exists (u°,c) G L^ (T) x 
[cq, Ci] which attains the supremum ( 13~T|) denoted by b: Let (uf, c l ) be the sequence for 
which || 0Aj( M ii cl ) — 0*(^i; c *) II z, 1 (T) — > b as i — > oo. Let i>° be a primitive of u° which 
belongs to Lzp r (T). We have a subsequence of (v®,c l ), still denoted by (w°,c 4 ), which 
converge to (v°,c). By Proposition 12.41 and \2.7\ we have Aj (u? x ; c l ) — >■ Aj (w°;c) i n 
L 1 (T) and 0*(u° x ; c l ) — > 0*(t>°; c) in L X (T) as i — > oo. Therefore we obtain || (j^.^x'i c ) ~ 

W;c) IUht)=&. 
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Let (u®,c>) attain (13.11) . Let v® be a primitive of which belongs to Lip r (T). We 
have a subsequence of (^,c J ), still denoted by (u°,o J ), which converges to (v°,c). It 
follows from 7. of Theorem 12.61 that there exists 5o > such that, if Ax < 5q, we have 

II 0aK; c ) -0*O°; c ) IIl i (t)< f • Hence 

y > II 0a 0°; c) - c) || L i (T) 

> II c*) - #(^; c*) |Ui( T ) - || A (^°; c) - A (^ ,; c*) || il(T) 

-|| 0*(^;c J, )-0 i K o ;c) || L1(T) . 

ByPropositionOandEZl || A O°; c)-0 A (^; C*) || L i (T ) + || c^)-0*(?;°; c) || L i (T) < 

Si for large j. Therefore we have || 4> A (Vj X ; c J ) — (^{y^x, cP) ||li(t)< £o f° r an Y Ax < 5 , 
which is a contradiction. □ 

Next we see that the convergence || (fi A (u°; c) — c) Hl^t) - ^ as A — )• 0, which is 

uniform with respect to (u°,c), yields uniform boundedness of || 4> A (u°;c) ||x,oo with the 
aid of the one-sided Lipschitz condition. 

Proposition 3.3. Let Ai > be that of Theorem \2.b\ with T = 1 and A = (Ax, At) be 
such that < Ao < A = At/Ax < Ai satisfying the conditions in Proposition \2. 8\ Then 
there exists 5 > such that, if Ax < 8, we have with j3\ in Proposition ^. 1\ 

SUp || 4>\{VP)C) ||i=o< + 1. 

«°6L^ (T),c6[co, Cl ] 

Proof. Let < e < 1 be such that ±=^5 > 2e T > / H* pp > E 2 A K , where 2 A' At = I. We have 
5 > in Lemma [3.21 with this e and t = 1. We take Ax < min{5, y/e}. Consider 

A := {y G T | |0 A (n°; c)(y) - ^(u ; c)(y)\ > y/E}. 

Since || A (w°;c) — 1 (m°;c) ||li(t)< £, meas [A] < y/e. Hence foryeA there exists 
x, x G K \ A such that < y — x < y/e, < x — y < y/e. For x G T \ A, we have 
|0 A (u°; c)(x) - ^(n ; c)(x)| < and |0 A « c)(x)| < \<P\u°; c)(x)| + ^(t) + 1. 

Consider 

A:={yeA \ \ ( j> 1 A (u°;c)(y)\>p 1 (t) + l}. 

Suppose that A is not empty. Then there exists x n G A n (AxZ) such that (i) > 
/?i(l) + 1 or (ii) it^ < -/?i(l) - 1. Since there exist x m ,x m / G (R \ A) n (AxZ) such 
that < x n — x. m < + 2Ax < 3 y/e, < x m / — x n < y/e + 2Ax < 3y/e, we have 



^ ~ < > /3i(l) + l-(/?i(l) + y^) = We > ^ 

3v^ 



, n - uj K -(/3 1 (l) + v ^)-(-/3 1 (l)-l) p2K 

— e Oyf e 

These two inequalities contradict the one-sided Lipschitz condition. □ 

Remark. 7. of Theorem 12.61 states that A (w°;c) converges to 1 (m°;c) uniformly on 
T \ as A —7- 0, where G is an arbitrary "small" neighborhood of shocks. We may 
not use this fact for Proposition 13. 3^ because it is not easy to obtain uniformness of the 
convergence with respect to (u°,c). 
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Proof of Theorem 13.11 Let Ai > be that of Theorem 12.61 with T = 1 and r > 
/?i(l)+l. Let 5 > be that of Proposition E3J Then, for any u° € L^ (T),u° := A (w o ;c) 
belongs to L^^^ (T). Hence, by the choice of Ai, A c) = c) is well-defined 

and bounded by /?i(l) + 1 again. In this way (f)\{u°] c) can be defined for t — > oo with 
the CFL-condition. □ 



3.2 Long Time Behavior 

If we take r > + 1, (Pa(u°; c) belongs to L^° (T). Therefore A ('j c ) * s a ma P from 
L£° (T) into itself. We can find fixed points of the map for each c. In this subsection, we 
consider the fixed points and their stability, which makes clear the long time behavior 
of the Lax-Friedrichs scheme. Note that the Lax-Friedrichs scheme has a contraction 
property under the CFL-condition: For < t < t' we have 

II 0*a(m°; c) - A (w°; c) |Ui ( t) < || 0a(«°; c) - a (m°; c) || l i (t) . 

This can be sharpened to be a strict contraction property. Let Ylm-k De summation with 
respect to {m | < m < 2N, m + k = even} for each fixed k and || x \\i'= Xa<j<n l x il 
for x £ W 1 . 

Proposition 3.4. The family of maps {0 A (-; c )}t>o has a strict contraction property 
within the time period 1: For any two different initial datas u° and u°, we have 

|| 0'+ V; c) - </4 + V; c) \W m <\\ c) - A («°; c) || L i (T) . 

Proof. It is enough to show that any two difference solutions u k ,u k of ( 12. 31) satisfy for 
all fc > 

|| u k+2K - u k+2K \\i<\\ « fc -w fc ||i. 
Set := — u^, cr^ := sign = 1 or — 1 (signO := 1). Then || u k — u k \\ x = 
J2m-,k \ z m\ — J2m;k a m z m- By tne difference equation of (12.31) . we have 

El ~k+l I _ V 1 _ ^k+l f*--k ( i \ rfc \ I 1 Jfc /-i i \ rfc \ 1 

l Z m+ll — / > a m+l Z m+l ~ / 4 a m+l| p m + 2 ^ A0 m+2J + 2 ~ l ~ AO m) } 1 

m\k m\k m;k 

where := H p (x m , tk, c + u k m + with a constant 6*^ derived from Taylor's formula. 
Shifting the order of the above sum, we obtain 

m;k m;k ^ 

= \ Z m\ l^rr; 

Let denote the second sum in the second line of the above equality. We see that 
R k < 0, since each term of [ ] in R k satisfies 

(1) If ^t\ + ^+\ = ±2,then[ ] = -l±^ = 0or-2, 

(2) If a k +\ + = 0, then [ ] = -1 ± < 0, due to the CFL-condition. 



-l + o- 



2 



4- -rr fc+1 I 
^ 2 m+1 



1 + A5* 
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Since u k ,u k have the same average equal to and u° ^ u°, z^ necessarily change the 
sign and the case (2) happens. It seems possible that even if u k ,u k are as such, we have 
R k = 0, namely z^ = for all the integers m for which the case (2) happens. But 
after some more fc*-time evolution (A;* < N < 2K) the case (2) certainly happens and 
R k+k * < 0, because such zero-points disappear as k increases due to the monotonicity of 
the Lax-Friedrichs scheme under the CFL-condition (see also Remark 3.4 in [13]). □ 

We now show the existence of time periodic difference solutions and their stability, 
which provides long time behavior of the Lax-Friedrichs scheme. 

Theorem 3.5. Take r > /?i(l) + 1 and fix A = (Ax, At) so that Theorem \2.6i and 
IJ.il holds. Then, for each c, there exists a fixed point u A G L^ (T) of (f) A (-;c), which 
yields a time periodic difference solution (p A (u A ;c). Such a periodic solution is unique 
with respect to c. Any other solution (p A (u°; c) exponentially falls into the periodic state, 
namely there exists p G (0, 1) and a > depending on A, but independent of u°, such 
that || (p A {u°; c) - <f) A {u A ; c ) IU°°< a P l f or t eN. 

Proof. The map 4> A (-;c) is actually a map from ~R N to M. N , since the step functions 
have only N different values at most. Let B r be the set of all x G M. N with || x ||oo< t. 
If r > /?i (1) + 1, then the map 4> A {-\ c) is actually a map from B r to B r . Therefore we 
obtain a fixed point due to a fixed point theorem. By Proposition 13.41 periodic solutions 
must be unique. Exponential decay can be proved in the same way as the proof of (5) 
of Theorem 2.1 in [T3]. □ 

Remark. It is likely in general that p gets arbitrary close to 1 as A tends to 0. 

Theorem 3.6. Take r > /?i(l) + 1 and fix A = (Ax, At) so that Theorem\EM andUl\ 
holds. Then, for each c, there exists a constant h&(c) G 1R such that if h(c) = h^(c) 
we have a fixed point G Lip r (T) of c) , which yields a time periodic difference 
solution V'aI^a! c )- Such a periodic solution is unique with respect to c up to constants. 
Any other solution ^ A (w°;c) exponentially falls into a periodic state, namely for the 
same p G (0,1) and a > in Theorem \3.5\ and d G E depending on (v°;c) we have 
|| ^(v °; c) — ip f A (v% + dh; c) ||c°< a P l f or t G N, where h(x) := 1 and ^(v + dl\\ c) = 
^ A (v°;c) + dh. 

Proof. We imitate the proof of the weak KAM theorem [8]. The map iPa(''> c ) ^ s actually 
a map from W N to R , since the piecewise linear functions have only iV different slopes 
at most. Let us write x ~ y for x, y G M N , if there exists b G M such that y = x + bl 
with 1 := (1, ■ • • , 1) G M. . Introduce x :— {y G M. N \ y ~ x}, || x \\:— inf^g^ || y ||oo, 
B r := B r j ~ and 

i)\{w]c) := y, 

where w,y G M. N are the vectors given as Wi := ip A (v°; c)(x 2 j_i), yi := ip A (v°; c)(x 2 j_i), 
% = 1, 2, 3, • • • ,N. Note that c) is actually a map from B r into B r . Hence we have 
a fixed point w satisfying ip A (w;c) = w. Therefore, for the linear interpolation of an 
element w G w denoted by v A , we have a constant b(c) G M such that 

v A = ^ A (v° A ;c) + b(c)h. 
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This relation means that v A yields a time periodic solution of ( 12. 5 j) with h(c) + 6(c) 
instead of h(c). 



Note that ^(w°;c) < ip A (v°;c), if t»° < Let a°,b° be constants such that for all 

x e T 



with at least one point attaining the equality in each inequality. Then we have v A (x) + 
b° < ip A {v°\ c)(x) < v A (x) + a for all x £ T. Let a 1 ,^ 1 be constants such that for all 



with at least one point attaining the equality in each inequality. Note that a 1 < a 
and b 1 > b°. Then we have v A (x) + b 1 < tp A (v°; c)(x) < v A (x) + a 1 for all x £ T. 
In this way we obtain the monotone bounded sequences a? and fr 7 . Take d such that 
Hindoo b 3 < d < lim^oo a? . Then v A + dii and ip A (v°; c) coincide at least one point for 
any t £ N. Let Xo £ T be such that v A (xo) + d = ifj A (v°; c)(xq). Then we obtain for all 
x £ T and t £ N 



Introduce the map h A {c) : c i— > /i(c) + b(c), which is the effective Hamiltonian of the 
difference Hamilton- Jacobi equation ( 11.61) . We remark that h A (c) may play an important 
role in numerical analysis of the weak KAM theory. It is meaningful to investigate its 
properties. 

3.3 Effective Hamiltonian 

Here is the characterization of h A (c), which is very similar to that of the effective Hamil- 
tonian h(c) of the exact Hamilton- Jacobi equations (ll.2p . We refer to [I] for the charac- 
terization of h(c). 

Theorem 3.7. 1. h(c) — h A (c) is the unique value for which M.b}) admits a space- 
time periodic difference solution. 

2. h A (c) is the averaged Hamiltonian: For the space-time periodic difference solution 
u k m of ( ti.5j) . we have 



v° A (x) + b° <v°(x) <v° A (x) + a { 







x £ T 



v° A (x) + b 1 < i[) A (v ; c)(x) < v° A (x) + a 1 




□ 



M c )= Yl J2H(x m ,t k ,c + u^(c)) -2AxAt. 



0<k<2K m;k 



3. Let ^ +1 (c) be a time global solution of the difference equation 
(3.2) D t v k + l + H{x m , t k , c + D x v k m+l ) = 0. 



Then we have for all n 




h A (c). 
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4- h/±(c) is a convex C l -function. 

5. h A (c) uniformly converges to the exact effective Hamiltonian h(c) of hi. 2^1 as A — > 
0: 

sup \Jia{c) — h(c)\ < f3 3 \/Ax. 

ce[co,ci] 



Proof. 1. Let be another space-time periodic solution of ( II. 6p with h(c) = h A (c). 
Extending periodic solutions to the entire odd-grid, we have the stochastic and varia- 
tional representation formula up to any negative time index Iq\ 



l <k<l+l 



V 



m(-y l ) 



E 



E ^(7 fc ,i fe -i,lt( 7 *))A^ 



m(7 i 0) 



+ h A (c)(ti +1 -t k ), 
+ ^a(c)(*i+i - ti„). 



Zo<A:<i+l 

By the variational property, we have 

=-i+i -i+i 



(3.3) < +i - < +i < £ M , r) 



m(7'o) 771(7*0) 



+ (/iA(c) - /iA(c))(t«+l - t k ). 



Noting that v^ +1 , v^ +l are periodic and hence bounded, dividing ( 13. 3ft by tj+i — 1; and 
letting Z —7- —00, we obtain 

< %a(c) - h A (c). 

A similar reasoning yields the converse inequality and concludes h A (c) = h A (c). 

2. Since v^+i satisfies DtV 1 ^ 1 + H(x m , t k , c + D x v^ +1 ) = Jia(c), we have 

Mc)= E Y. D ^™ 1 - 2AxAt+ E E^^-'^^+^+i)" 2 ^^- 

0<k<2K m;k 0<k<2K m;k 

The first term of the right hand side is equal to due to periodicity of v^ +1 . 

3. Let t^ +1 (c) be the solution of DtV^ + Hfam, tk, c+_D x {;^ +1 ) = h A (c) with the same 
mesh size as ( 13. 2 j) and = w m+i- ^ follows from Theorem 13.61 that we have, adding 
a constant if necessary, |{^ +1 (c) — ^ +1 (c)| — > as Z — > 00. Therefore is bounded for 
I — > 00. Since 



inf E 



E 

'0<k<l+l 



L c ( 7 fc ,t fc „ 1 ,e (7fc) )At + ^ (7 o ) 



mfE M( , [ mi^h-i,e mi , k) w+vi { J+h A ( C )t l+1 



0<k<l+l 



and their minimizing velocity fields are the same, we obtain i^ +1 (c)— ?)^ +1 (c) = — /ia(c)£/+i. 

4. Following the proof of (6) of Theorem 2.1 in [13], we can prove that c + ^(c) is 
a C 1 -function for each m,k. Therefore 2. yields (^-regularity of h A . Let w^ +1 (c) be a 
solution of (13.21) . Fix n. The map c 1— > t^ +1 (c) is a concave function for each I + 1 > 1: 
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Let £* be the minimizing velocity field for f^ +1 (c*) with c* := Ac + (1 — A)c, A G [0,1]. 
Then we have 



^ +1 (c*)-{A^ +1 (c) + (l-AK +1 (£)} 



0<fc<Z+l 



;i-A)^ ( . ;e) [ £ -(c*-c)f 



fc Ay- 



0</c<Z+l 



0. 



Therefore the map c H- w„ (c)/t/ + i is also a concave function and h A (c) = — lim^oo v + (c 
is a convex function. 

5. Hereafter &i, &2, • • ■ are positive constants independent of A and c. For each x G T, 
we have m such that x G [x m+ i,x m+3 ). Note that v^ +1 < v A (x, 1) < w^? fl with 
(n*, n*) = (m, m + 2) or (m + 2, m) and 

(3.4) _ v(x nt+1 , 1) - 2rAx < v A {x, 1) - 1) < v 2 J +1 - v{x n . +1 , 1) + 2rAx. 

Let x G T attain max 2/eT (vA(?/, 1) — v(y, 1)) and n* be defined in the above way with 
the x. Let 7* be a minimizing curve for v(x n * + i, t). Define £ as £^ := 7*'(£fe). Note that 
?7 fc (7) defined by the £ satisfies 1 77^(7) — 7*(^)| < b\ Ax for any < k < 2K. By the 
representation formulas and Proposition 12. 5[ we have 



v(x n . +1 ,l)= / L c ( 7 *( S ), S ,7* / ( S )) C / S + t;(7*(0),0) + /i(c), 
./o 

(3.5) ^f +1 <E K , [ £ L c (7 fc ,t fc - 1 ,e (7 M) At + ^(7°,0)] + fc A (c) 

o<fc<2_ft: 

]T L c (r / fc (7),t fc - 1 ,^ (7 , ) )At + t; A (r / (7),0)] + ^ A (c) + 6 2 VA^ 

<fc<2_ft: 

E L c (7*(4),t fc _ 1 ,7* / (t fc ))At + ^ A ( 7 *(0),0)j +h A (c) + b 3 yfAx~ 



< E, 



M(-iO 



L 0<fc<2K 

< / L c ( 7 *(s),s,7* / (s))c/s + t; A (7*(0),0) + / iA (c) + 64V / A^. 

Therefore, noting (I3.4j) . we have 

w A (x, 1) - u(x, 1) < u A (7*(0). 0) - u(7*(0), 0) + h A (c) - fc(c) + b 5 VAx~. 

It follows from the periodicity of v A , v and the choice of the x that (v A (x, 1) — v(x, 1)) — 
(«a(7*(0),0) - w(7*(0),0)) > 0. Therefore we obtain 

—b 5 V~A~x < h A (c) — h(c). 

Let x G T attain min ?/g x(u A (y, 1) — ?>(?/, 1)) and be defined in the above way with 
the x. Let be the minimizing velocity field for v 2K +1 . Then we have 



v 



2K 
n.+l 



E^J ^(7 fc ,4-i,et (7 ,))At + ^ A (7 ,0)j + h A (c), 



0<k<2K 



> £ L c (r / fe ( 7 ),4-i,et (7fc) )At + t; A (r / (7),0) 



+ ^a(c) — fefivAr. 



0<k<2K 
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Let 77a(t) be the linear interpolation of r] k (j). Note that ^a(7)'(0 = ^*m( 7 fe ) ^ or * e 
(tk-i,tk). For each 7 we have 

(3.6) ^Ow, 1) < / L%7 ]A ( 1 )(s),s, VA ( 1 )\s))ds + v( VA ( 1 )(0),0) + h(c) 

Jo 

< Yl iC ^W>Vi, r^))At + tXr/°( 7 ), 0) + h(c) + 6 7 Ax. 

0<k<2K 

Therefore, noting (13. 4p . we have 

v A (x, 1) - v(x, 1) > £ M , e) [^(^(t), 0) - v(r}°(j), 0)] + h A {c) - h(c) - b 8 VAx~. 
Since (v A (x, 1) — u(x, 1)) — (^a^ ^)? 0) — v^rj ^), 0)) < for all 7, we obtain 

^a(c) - h(c) < b 8 VAx. 

□ 

3.4 Convergence of Periodic Solution 

It is easy to prove convergence of space-time periodic difference solutions to exact ones 
up to subsequence. We should remark that space-time periodic viscosity solutions and 
entropy solutions are not unique with respect to c in general. The selection problem in 
finite difference approximation remains open. It is also challenging to investigate details 
of the convergence even in the case where the uniqueness holds. We will make some 
progress in this issue in the next section. 

Theorem 3.8. There is a sequence A = (Ax, At) — > for which {v A } and {u c A } 
converge to a Z 2 -periodic viscosity solution v of ( dUP with h(c) = h(c) and a Z 2 -periodic 
entropy solution u = v x of U.l\) respectively: 

SUp || V° A (-, t) - «(-, t) \\ C a-> 0, SUp || U° A (-, t) - U(-, t) ||£i(T)-» 0. 

ieT teT 

Proof. Since {v A (-, 0)} is a family of Lipschitz functions with a common Lipschitz con- 
stant, we have convergent subsequence, still denoted by v A (-, 0): v A (-, 0) — > v°. Let v be 
the viscosity solution of (II. 4p with v° = v° and h(c) = h(c). Then we have a minimizing 
curve such that 

v{x n , t l+1 ) = [ L c ( 7 *(s), s, y*'(s))d8 + ^°(7*(0)) + h(c)t l+1 . 
Jo 

By an estimate similar to f)3.5p . we have 

v A (x n , t l+1 ) < [ L c (7*(s), s, y*'{s))ds + v A (y*{0), 0) + h A (c)t l+1 + hVAx'. 



o 



Since h A (c) — > h(c), we obtain 



limsup{u A (x n ,t /+ i) -v(x n ,ti +1 )} < 0, 

A->0 
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which is uniform with respect to {x n ,ti +1 ) G T 2 . By an estimate similar to (13. 6p . we 
obtain 

liminf{w A (x n ,t m ) - v(x n ,t l+1 )} > 0, 

A— >0 

which is uniform with respect to (x n ,ti+i) G T 2 . Therefore we conclude that v A — > v 
uniformly on T 2 and v is Z 2 -periodic due to periodicity of v&. 

It follows from reasoning similar to the proof of Theorem 2.8 in [T8] that u c A : = 
(^a)z converges to u = v x pointwise almost everywhere in T 2 , where {v A } is the above 
convergent subsequence. Hence we have || u%(-, t) —u(-,t) > for each t. It follows 

from reasoning similar to the proof of Proposition 2.14 in [13] that u c A satisfies 

II w A (-, tk) — u c A (-,tk') ||li(t)< b 2 \tk — tk<\ 

with a constant b 2 independent of k, k', c and A. Therefore u G Lip(T; L 1 (T)) with the 
Lipschitz constant b 2 . Thus we conclude the theorem. □ 

4 Error Estimate 

We show error estimates for entropy solutions of initial value problems and for Z 2 - 
periodic entropy solutions in the special case where they are associated with KAM tori. 
The latter is a rigorous result on finite difference approximation of KAM tori. We refer 
to [2] for an error estimate for Z 2 -periodic entropy solutions associated with KAM tori 
in the vanishing viscosity method. 

4.1 Error Estimate for initial value problem 

The following error estimates hold: 

Theorem 4.1. Let A = (Ax, At) satisfy the conditions in Theorem \2.6i and Proposition 
\2.8i . Let u be the entropy solution of ( fj.^j) and ua be given by the difference solution of 
A2.3\) . Then the following holds: 

1. For any T G (0, oo) and each t G (0, T\, there exists a constant fi±(t) > indepen- 
dent of initial datas for which we have 

II u A (-,t) - u(-,t) ||li(t)< /3 4 (t)Aa:z. 

In particular, ifu° satisfies the one-sided Lipschitz condition u°(x)—u°(y) < M(x— 
y) for a.e. x > y and M > 0, then there exists a constant f3 5 > for which we 
have 

SUp || U A (-,t) - u(-,t) 1 1 Li (T)< A) Ax 
0<t<T 

2. If u is Lipschitz mix [0, T], then there exists a constant f3 6 > for which we have 

sup |ma(x, t) — u(x, t) | < /3 6 Ax * 

O,t)6Tx[0,T] 
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Proof. 1. By Theorem 12 .6[ we have for all t G [0, T] and all initial datas 

II v A (-,t) - v(-,t) \\ C o< f3 2 VAx. 



By Proposition I2.8[ we have for each t G [At, T] and all initial datas 



(4.1) 



U A (x m+2 , t) - U A (x m , t) 

2Ax 



< E 



k(t) 



Since u A (-,t) has the average 0, we have 

— u A (x m , t)} = ^2{u A (x m+2 ,t) — u A (x m , 

t)} 

m;k(t) m:+ 

+ t) ~ U A (x m , t)} 

m: — 

= o, 

where J2 m -+ (Em:-) stands for the summation with respect to m for which u A (x m+ 2, t) — 
u A (x m ,t) > (u A (x m+ 2,t) — u A (x m ,t) < 0). Hence it follows from (14. ip that the total 
variation of u A (-,t) on T is bounded: 



(4.2) ^2 \ u A(x m+ 2,t) - u A (x m ,t)\ = 2y^ y {u A (x m+2 ,t) - u A (x m ,t)} < 2E 

m;k(t) 



k(t) 
A • 



m: + 



For any e > 0, there exists A = (Ax, At) such that 

II u A (-,t) - u(-,t) ||li( T )< e. 

In particular, we take such a A = (Ax, At) that satisfies At/ Ax = At/ Ax, Ax < 
(/^e) 4 , Ax/ Ax = 3 P for some p G N. The last relation guarantees that the points of 
discontinuity of ha are also these of u A . Then we have 



(4.3) 



U A (-,t) - u(-,t) \\ L i {T) < || U A (-,t) - U A (-,t) +£, 



t)-v A {-,t)\\ c o < / 3 2 v / A^+ftVAi<2ftv / A^. 



Now we estimate || u A (-,t) — u A (-,t) \\l x (t)- We introduce w A := u A (-,t) — u A (-,t), 
w A := v A (-,t) — v A (-,t) and k(t) := 3 p k(t). Let x m G AxZ and x mo := for k(t) =even 
Ax for k(t) =odd. We divide AxZ according to the sign of w A : Define 



or x 



/Do 



i In+i as 

h 
h 
h 

In 
In+1 



{x mo ,x mo+ 2, ■■■ ,x mi } on which w A (x) > (< 0), 
{x mi+ 2,x mi+4 , ■ ■ ■ ,x m2 } on which w A (x) < (> 0), 
{x m2+2 ,x m2+4 , ■ • • ,x ma } on which w A (x) > (< 0), 

{^m n _ 1 +2,^m n _ 1 +4, ' " ' , X m „ } Oli which W A (x) < (> 0), 

{x mn +2, x mn+i , ■■■ , x mo + 1} on which w A (x) > (< 0), 



where n =even and x m „ < x 



2Ax. We redefine Ji as 



I\ . — {x mo ,X mQ +25 ' ' ' i^mj With X mQ . — X rrin J r 2 1. 
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Note that w A (x) > (< 0) on Set := x mi — x mo + 2Ax and := x mj — 
^mj_i+2 + 2Ax for j > 1. We have Y^=i l-^il = 1- For each J 3 - on which > (< 0), 

we have yi G 7j for which w A (x) takes the maximum (minimum) within Ij. Suppose 
that w A (x) > on I\. In the other case, the argument is parallel. Note that 

n/2 ( 

II w A (x) \\li(t)= ^ < w)a(i) • 2Ax - ^ «)a(i) • 2Ax 

j=i ^xei"2j-i a;e/2j 

Introduce J := {j | < j < n/2, max{|/ 2i _i|, \I 2j \} < Ax 1 / 4 } and J := {j | < j < 
n/2, max{|J 2i _i|, \I 2j \} > Ax 1/4 }. We have jjj- Ax 1/4 < 1 and jjj < AoT 1 / 4 . Therefore, 
noting w;a = (w A ) x , (14. 2[) and (14. 3p . we obtain 

II W A (x) ||ii(T) = 1 ' ~~ ' 

+ ' 2 ^ x - X] w a(x) • 2Ax 

< Yl \ w ^(y 2j ^) ~ w A (y 2j )\Ax^ 

+ ^ {^(2^.1 + Ax) - W A (x m2j _ 2 +2 ~ Ax)} 

-{w A (x m2j + Ax) - l5 A (Vi+2 ~ 

< (2x^ (t) + 2x^ {<) )A^ +ttJ-4-2/3 2 v / Ax~ 

< 4E^ w Ax3 + 8/3 2 Ax^, 
Since e is arbitrary, we conclude that 

II u A (;t)-u(-,t) \\ L i m < {AEf ] + 8/3 2 )Axi 

If t>° satisfies the one-sided Lipschitz condition, < max{M, E*} for all t. 

2. Note that it = v x and u is C 1 ' 1 with |t> A (x, t) - v(x, t)\ < (3 2 y/Ax on T x [0,T]. Fix 
arbitrary i 6 [0, T]. We have for any x, x' G T 

«a(2/, <) - < 2/3 2 v / Ar. 



We already know that the step function u A (-,t) satisfies the one-sided Lipschitz con- 
dition. Therefore the linear interpolation of Um with respect to the space variable, 
denoted by u A (x), satisfies || u A — u A (-, t) ||z,i(t)< feiAx. Setting w A := u A — u(-, t), we 
have for all x, x 1 G T 

(4.4) | f w A (y)dy\ < b 2 ^SAx~. 

Jx' 

Since u is Lipschitz, w A still satisfies the one-sided Lipschitz condition 

w A (xi) - w A (x 2 ) 



< 6 3 . 



x 1 - x 2 
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Note that wa does not necessarily satisfies the Lipschitz condition, because «a does 
not necessarily satisfy the Lipschitz condition. Suppose that |ioa(^)| > b^Ax* with 
(6 4 ) 2 /(46 2 ) > ^3 for some x. Let / 3 x be a connected interval on whose boundary we 
have |u>a(^)| = |Aii. By (I4.4p . we see that 

|J| < ^Axi 

If w A (x) > (< 0), we have with the left (right) boundary of / denoted by x 

w A (x) - w A (x) > (6 4 ) 2 { w A (x) - w A (x) > (6 4 ) 2 

x — x ~ Ab 2 3 V x — x ~ Ab 2 

which is a contradiction. Therefore we obtain 

II W A \\c°< ^Ax*. 



Since \u A (x,t) —u(x,t) \ = — u(x, t)\ < \um^ —u(x m ,t)\+b 5 Ax = \w A (x m )\ +b 5 Ax, 
we conclude the theorem. □ 



4.2 Error Estimate for KAM Tori 

Let u c = v c x be a Z 2 -periodic entropy solution of the C 1 -class. We remark the relation 
between such a u c and Hamiltonian dynamics: Consider the time-1 map f : TxM. — > TxM. 
of the Hamiltonian flow generated by the flux function H(x,t,p) with the initial time 
equal to 0. Then {(x, c + u c (x, 0)) | x G T} = T is a smooth invariant torus of /. Due 
to the classical result of Poincare, there exists a rotation number oj\. Let us regard the 
non-autonomous Hamiltonian dynamics generated by H(x,t,p) as the autonomous one 
generated by 'H(qi,q2,Pi,P2) '■= P2 + H(qi,q2,Pi) in the extended phase space T 2 x M 2 . 
Define 

l(u c ):={(q,g(q))\q=(q u q 2 )eT 2 }^T 2 , 

where g(q) := (c + u c (qi,q 2 ),h(c) — H(q 1 ,q 2 ,c + u c (q 1 ,q 2 )). Then X{u c ) is a smooth 
invariant torus of the Hamiltonian flow tpifa generated by Ti. Let C(s) := (7*(s),s) 
be the characteristic curves of u c , which satisfy 7*'( s ) = -^p(7*( s )) s,c + u c (7*(s), s)) for 
sGl. The dynamics of the reduced characteristic curves C*(s) := C(s) mod 1 = (7*(s) 
mod 1, s mod 1) and that of the trajectories on X(u c ) are identical, namely we have for 
all s e R 

cp s n (C*(0),g(C*(0))) = (C*(s),g(C*(s))). 

Due to the classical result of Poincare, C(s)/s converges to cu — (cji, 1) G M 2 independent 
of C(0) as |s| — > 00. This 00 is called a rotation vector of T(u c ). If the rotation vector 
is irrational, each trajectory starting from a point of X{u c ) is dense on I(u c ). Therefore 
we obtain information on u c from only one characteristic curve, which is the key idea of 
the following argument. Approximation of u c leads to that of the invariant torus I(u c ). 

Now we consider a special case where X(-u c ) is a KAM torus: We say that c is associated 
with a KAM torus, if u c is C 1 and the dynamics of C* (s) is C 1 -conjugate to that of a linear 
flow on T 2 with a Diophantine rotation vector, namely there exists a diffeomorphism 
F : T 2 -> T 2 such that 

C*(s) = F(us + 6), 
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where Sel depends on C* (0) and u G R 2 satisfies the v, r-Diophantine condition 
\u\Z\ + > v || z \[ T f° r a H z £ Z 2 \ {0}. 

If c is associated with a KAM torus, u c is the unique Z 2 -periodic entropy solution of 
( 11. ip with the c. For the existence of a value c associated with a KAM torus, we 
refer to the KAM theory to the autonomous Hamiltonian systems with two-degree of 
freedom generated by the above H(qi, ?2>f>i>£>2)- We remark that the KAM theory under 
Russmann's non-degenerate condition (see e.g. [16]) works for such a degenerate % in 
(p!,p 2 ) with additional assumptions. The following error estimates hold: 

Theorem 4.2. Let A = (Ax, At) satisfy the conditions in Theorem \2.6] Proposition 
\2.8i and Theorem \3.1\ Suppose that c is associated with a KAM torus. Let v c be a 
T? -periodic viscosity solution such that v c x = u c . Then we have the space-time periodic 
difference solutions v c A and u c A that satisfy 

sup \v A (x,t) - v c (x,t)\ < /3 7 Ax^+^, sup \u c A (x,i) - u c (x,t)\ < (3 8 Ax*^ , 
(x,t)eT 2 (x,f)ex 2 

where j3 7 and /?§ are independent of A. 

Proof. Let v c A be a periodic difference solution. We omit c in v c , u c , etc. Fix an arbitrary 
t 6 T. By adding a constant to v&, if necessary, it holds that v&(x,t) — v(x,t) < for 
all x G T and va(x*, t) — v(x*,t) = for some x* G T. Then we have n* and / such that 

= i'A{x*,t) - v(x*,t) < v l n » +1 - v(x n *+i, U) + 61 Ax, 

where \x* — x n *+±\ < 2 Ax and t G [ti,ti + i). For any j G N, we have a minimizing curve 
7* such that 

v(x n . +1 ,ti) = [ L c { 1 *{s),s, 1 *'{s))ds + v{ 1 *{-j +ti),-j + t l ) + h{c)j, 
J-i+ti 

v l n * +1 < / ^ C (7*(s), s, l*'(s))ds + v A (>y*(-j + t t ), -j + U) + h A (c)j + b 2 VAxj, 

where we use an estimate similar to ( 13. 5p . Hence we obtain with 5. of Thorem 13.71 

< v l n , +1 - v(x n * +1 ,ti) + hAx 

< v A (l*(-j + tt), -j + ti) - v(l*(-3 + U), -J + tt) 
+ {h A (c) - h(c))j + b 2 VAxj + bi Ax 

< v A {l*H +t),t)- vtfi-j +t),t) + b 3 VAx~j. 

Since v A (x, t) — v(x, t) < for all x G T, we obtain 

\v A (i*(-j + t),t) - v(j*(-j + t),t)\ < hVtej. 

Since C*(— s + t) := (7*(— s + t), —s + t) mod 1 is a reduced characteristic curve, we 
have C*(-s + 1) = F(w(-s + t) + 9). It follows from [6] and |3] that for the set 

A4 := {9 + co{-s + 1) mod 1 I < s < ^} 
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is e-dense on T 2 , namely 

IJ B e(0 = T 2 , 

where B £ (() = {( G T 2 ||C - C ||i< e}- Define T := F^fT x {*}) and X := (x,t) for 
x G T. For each X, we have ( G JV e fl T such that || X - C ||i< £ with X := F~ 1 (X) 
and ( = co(— s* + t) + 6 mod 1 with some < s* < Note that s* must be an integer, 
because F(() = C*(—s* + 1) G T x {t} and — s* + t mod 1 — t. Hence, denoting s* = j, 
we have 

II X - C*(-j + t) HHI *W " ^(0 lli<ll DF \\ op e. 
Therefore we obtain for all x G T 

\v A (x,t)-v{x,t)\ < \v A (F(X))-v A (F(())\ + \v A (F(())-v(F(0)\ 

+\v(F(Q)-v{F(X))\ 

< b 5 e + b^y/Axj + b 5 e 
^ . ,v / Ax . 

< h(— — + £). 

1 

Taking e = Ax 2 < 1+r ) , we have for all x G T 



\v A (x, t) — v(x, t) I < 2o 6 Ax 



2(l + r) 



Note that 66 is independent of the choice of t. We have for u A = (v A ) x , u = v x and all 
x, x' G T 



u A {y, t) - u{y, t)dy\ < 4o 6 Ax 2 d+-> . 

Since u A satisfies the one-sided Lipschitz condition, the linear interpolation of u l m with 
respect to the space variable, denoted by u A (x), satisfies || ma — u A (-,t) ||l1(t)< b^Ax. 
Setting w A := ma —u(-,t), we have for all x, x' G T 

(4.5) | / w A (y)dy\ < OsAx^ 1 ^. 

Jx' 

Since u is C 1 , w A still satisfies the one-sided Lipschitz condition 

w A (xi) - w A (x 2 ) 

S Og. 

x 1 - x 2 

Suppose that |wa(^)| > OioAx 4 ( 1+T ) with (6 10 ) 2 /(46 8 ) > b 9 for some x. Let / 3 x be a 
connected interval on whose boundary |iy^(a;)| = ^ a Aa; 4 ( 1 + T ) . By (14.5 j) . we see that 

|/| < — -Ax^+t). 
bio 

If wa(S) > (< 0), we have with the left (right) boundary of / denoted by x 
w A (x) - w A (x) (b w ) 2 ( w A (x) - w A (x) (6io) 2 , 

■> — — > Oq I -> — — > Ot 



x — x 46 8 V x — x 46 8 

which is a contradiction. Therefore we obtain 

II w A \\ c o< hoAx*^. 

Since \u A (x, t) - u(x, t) | = \u l m - u(x, t) | < \u l m - u(x m , t) | + b n Ax = \w A (x m ) | + b n Ax, 
we conclude the theorem. □ 
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